Evolution of Primordial Black Holes in Loop Quantum Gravity 
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In this work, we study the evolution of Primordial Black Holes within the context of Loop Quan- 
tum Gravity. First we calculate the scale factor and energy density of the universe for different 
cosmic era and then taking these as inputs we study evolution of primordial black holes. From our 
estimation it is found that accretion of radiation does not affect evolution of primordial black holes 
in loop quantum gravity even though a larger number of primordial black holes may form in early 
universe in comparison with Einstein's or scalar-tensor theories. 
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. I. INTRODUCTION 

X 

d ' The demand for consistency between a quantum description of matter and a geometric description of space-time 
indicate the necessity of a complete theory of quantum gravity. This theory is expected to provide a new light on 
singularities present in classical cosmology. Einstein's theory of general relativity leads to the occurrence of space-time 
singularities in a generic way. So, one may say, General Relativity is severely incomplete and is unable to predict 
what will come out of a singularity. One of the outstanding problems in classical Einstein cosmology is the Big- 
Bang singularity which is expected to be solved by quantum gravity. Loop Quantum Gravity (LQG) Q is one of 
the best motivated theories of quantum gravity. LQG is a background independent, non perturbative approach to 
quantum gravity. When Loop Quantum Gravity is applied to cosmology to analyse our universe, it is called Loop 
Quantum Cosmology(LQC) Q (also see [1] for a comprehensive review on LQC) . In loop quantum cosmology, the 
non perturbative effects add a term of -p^/pc to the standard Friedmann equation Q, where p represents the 
energy density of the universe and pc is the critical density at which the universe is completely filled with a free 
massless scalar field when the scale factor reaches a minimum. The modification becomes important when the energy 
density of the universe approaches critical density(pc) and causes the quantum bounce. So the classical big bang is 
replaced by a quantum big bounce in such a quantum theory of gravity. Recently more and more researchers have 
paid their attention to LQC inspired by its appeaiingieatures, like avoidance of various singularities Q, inflation in 
loop quantum cosmology(LQC) f^, large scale effect |8|, present cosmic acceleration [9| and so on. 
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Primordial Black Holes(PBHs) are the black holes formed in the early universe [10|. A comparison of the cosmolog- 
ical density at any time after the Big Bang with the density associated with a black hole shows that PBHs would have 
mass of the order of the particle horizon mass at their formation epoch. Thus PBHs could span enermous mass range 
star ting from 10~^ g to the typical values of lO'^^ g. These black holes could be formed due to initial inhomogeneities 
[ill [T2f. inflation [ll,[l3l) phase transitions [l^, bubble collisions [H, [13 or the decay of cosmic loops [ll|. In 1974 
Hawking discovered that the black holes emit thermal radiation due to quantum effects So the black holes get 
evaporated depending upon their initial masses. Smaller the initial masses of the PBHs, quicker they evaporate. But 
the density of a black hole varies inversely with its mass. So high density is required to form lighter black holes and 
such high densities are available only in the early universe. So primordial black holes are the only black holes whose 
masses could be small enough to have evaporated by present time. There have been speculations that PBHs could 
act as seeds for structure formation [20] and could also form a significant component of dark matter [2l|. Since the 
cosmological enviornment was very hot and dense in the radiation dominated era, an appreciable amount of energy- 
matter from the surroundings can be absorbed by black holes. Such accretion is responsible for the prolongation of 
life time of PBHs (Hill. 

In this work, we study the evolution of PBH within the context of loop quantum gravity. First, we estimate the 
cosmic scale factor a(t) and energy density p (t) of the fluid filling the universe for different cosmic era within the 
context of loop quantum gravity. Taking these as inputs, PBH evolution is studied considering both the Hawking 
evaporation and accretion of radiation by the PBH. The primary aim being to compare the results so obtained with 
the analyses carried out earlier within the context of General Theory of Relativity and Brans-Dicke theories. 



II. SOLUTION OF FRIEDMANN EQUATIONS 

For a spatially flat FRW universe(k=0) with scale factor(a), two Friedmann equations in loop quantum gravity (sj, 
take the form : 



\a/ 3 \ Pc 

2p 



H= ^ P + P (l- 

3 ^ Pc 

where H is the Hubble parameter, p is the energy density and P is the pressure of the fluid filling the universe. 
The energy conservation equation is given by 



(1) 
(2) 



p + 3H{p + p) =0 

From energy conservation equation, we get 



(3) 



p oc a for t < te 
p oc a^^ for t > te 

where te is the time of radiation-matter equality. 

Using this solution in equation(l), one gets the temporal behaviour of the scale factor a(t) as shown below. 
For radiation dominated era {t < te) : 



(4) 
(5) 



a{t) 



SttG, 



Pc 



3„ s 1/2, 2t 1/4 



) }■ 



(6) 



where the subscript indicates the present value of any parameter and = a(te) and pc represents the critical 
value of energy density of the universe given by pc — ig^^s Ppi with 7 — is the dimensionless Barbero-Immirzi 
parameter [24] and ppi is the energy density of universe in Planck time. 

For matter dominated era [t>te) : 



a[t) 



Poao' 



Pc 



1/2-1 2 



, r3 1/2 3/2 /SttG , n , / 3 Po 3\''^\ 
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(7) 



Using equations (4) and (6), we get 
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For t <te 



Pit) = PO + {^\l^P0'^'it - te) + l^P0"\t^ - to) 

1/2|2-| -1 



Using equations (5) and (7), we get 
For t>tp 



Pc' 



(8) 



(9) 



III. ACCRETION OF RADIATION 

When a PBH evolves through radiation dominated era, it can also accrete radiation from the surrounding. The 
accretion of radiation leads to increase of its mass with the rate given by 



Mace = 4:T^frBH Pr 



(10) 



where Pr is the radiation energy density of the surrounding of the black hole, rsH is the black hole radius and / is 
the accretion efficiency. The value of the accretion efficiency / depends on the complex physical processes such as 
the mean free path of the particles comprising the radiation surrounding PBHs. Any peculiar velocity of the PBH 
with respect to the cosmic frame could increase the value of / [25l [26| . Since the precise value of / is unknown, it is 
customary [l^l to take the accretion rate to be proportional to the product of the surface area of the PBH and the 
energy density of radiation with / ^ 0(1). 

After substituting the expressions for tbh = 2GM and pr given by equation(8) in equation (10) , we get 



Mace = WirfGHl^po 



Solving equation(ll), one can find 



Pc ^ 



8nG 1/2 



p'„'-{t-tc) + 



SttG 1/2, , 

— — Po' [tc - to) + 



On J 



M(t) = {Mr^ - S^fG^Pc''^^f^tar^-^[^^Z{t) if/^ + S^/G^.i/^^ 

tan-^[^Z{U)-lY/^\ 
Po J 

Again using horizon mass as initial mass of PBH i.e; Mi ~ MH{ti)=G^^ti, we get 



8ttG 

-1 



(11) 



(12) 



M{t) = M,[l - 8nfG'/'pc'/hJ-f{tan~^[^Z{t) 1]^^' tan-^[^Z{t,) 1]^/'}]-' 

V Stt Po Po 



where Z{t) = + {2J^po'l\t - Q + y^Po"^{tc - to) + (1 - f)''^} 



(13) 



IV. EVAPORATION OF PBH 



As is well known black holes can also loose mass through Hawking evaporation. The rate at which the PBH mass 
decreases due to evaporation is given by 
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Mevap = -^-irrBHO-HTBH 



(14) 



where an ~ is the Blackbody Constant and Tbh ~ is the Hawking Temperature = g^^Q^ 
Now equation(14) becomes 



M, 



an 1 



2567r3 G2M2 



Integrating the above equation, we get 



M = 



+ 3a{ti-t) 



1/3 



where a = ^ ^ 

and Mi is the the black hole mass at its formation time U 
We rewrite equation (16) as 



(15) 



(16) 



M{t) = Mil + ^{U - t) 



1/3 



(17) 



V. PBH DYNAMICS IN DIFFERENT ERA 



Primordial Black Holes, as discussed earlier, are formed only in radiation dominated era. We now study PBHs so 
formed in two categories: 

(i) PBHs evaporating in radiation dominated era {tevap < te) 
(ii)PBHs evaporating in matter dominated era {tgyap > te) 

CASE-1 (terap < te) 

If we consider both evaporation and accretion simultaneously, then the rate at which primordial black hole mass 
changes is given by 



+ 



3 8ttG 



MpBH 

"\t, - 



-Po 



m-KfC^M^Po 

to) + (l - 
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SttG 



Po'/'(t - te) 



rr - 



1 



2567r3 G2M2 



(18) 



Solving above equation numerically, we construct the following table for PBHs which are evaporating in radiation 



dominated era. In our calculation, we have used pa = 1.1 x 10 g — cm 
6.673 X IQ-^dyne - cw?/g'^, U = 10"s, U = 4.42 x lO^^s and M^ = lO^^g. 



5.317 X lO^^g - cm-^ G 



u 




tevap 






f = 


/ = 1 






3.333 X 10-"s 


3.333 X 10-"s 






3.333 X 10-''s 


3.333 X 10--'s 






3.333 X 10^ s 


3.333 X lO^s 


10-26s 




3.333 X lO^s- 


3.333 X lO'^s 



TABLE I: Display of formation times and initial masses of the PBHs evaporating in radiation dominated era 



It is clear from Table-1 that with increase in initial mass, evaporating time increases. However radiation accretion, 
surprisingly, seems to have little effect on evolution of PBH unlike the results obtained in theories of Einstein or 
scalar-tensor type. This is also shown in Figure-1. 
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FIG. 1: Evaporation of PBHs for different initial masses (i.e. lO^g, l(fg, W°g and lO^^p) are shown in tlie Figure where axes 
are taken in logarithmic scale 

CASE-II [tevap^ > 

Since there is insignificant accretion of radiation in matter dominated era, the first term in the combined equa- 
tion(18) for variation of Mpbh with time should be integrated only upto tg. Basing upon the numerical solution we 
construct Table-2 for the PBHs evaporating by the present time. 



tevap = to= 4.42 X 10"s 


/ 


U 


Mi 





2.3669 X lO^^^s 


2.3669 X lO^^g 


0.25 


2.3669 X 10"^^s 


2.3669 X W^^g 


0.5 


2.3669 X IQ-^^'s 


2.3669 X lO^^g 


1.0 


2.3669 X 10"^^s 


2.3669 X W^^g 



TABLE 11: Display of formation times of PBHs which are evaporating now for several accretion efficiencies 
It is clear from Table-2 that PBH evaporation is again not affected by radiation accretion efficiency. 



VI. CONSTRAINTS ON PBH 

The fraction of the Universe mass going into PBHs at time t is [9] 



m = 



(1 + ^)- 



(19) 



where Q,pBH{t) is the density parameter associated with PBHs formed at time t, z is the redshift associated with 
time t and Via is the microwave background density having value 10"^. 
Substituting the value of £7^ in the above equation, we get 



For t < te, redshift definition implies 



m = {i + zy^npBH{t) X 10^ 

a{t) a{t) a(tf.) 



(1 + ^)-' 



a(to) a{te) a{tQ) 



(20) 



(21) 
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But here 



a{te) 



Pc al 



3/2 



3/2 



3 V pca-^J ) . 



1/4 



and 



a{te) 
a{to) 



Pr 1^ 2 



1/3 



= 0.746 



Using above numerical values of different quantities in equation(23), we get 

a (to) 

Using equations (22) and (24) in equation (21), we get 

3/2 

-Pcai L 

Substituting equation(25) in equation (20), we get 

,3 , „3/2 
-Pc a\ 



(1 + ^)-' = 



Pc a; 



1/4 



X 0.746 



/2 



SttG 



PcaiJ i 



1/4 



X 0.746 X npBH{t) X lO"* 



(22) 



(23) 



(24) 



(25) 



(26) 



Using M = G ^t, we can write equation (26) to represent the fraction of the Universe going into PBHs' as a 
function of mass M as 



3/2 



SttG 



G{M - Me) + 1 



Peal) i 



1/4 



X 0.746 X npBH{M) X 10" 



(27) 



Observations of the cosmological deceleration parameter imply that Qbh{M) < 1 over all mass ranges for which 
PBHs have not evaporated yet. But presently evaporating PBH s{M^,) generate a 7- ray background whose most of 
the energy is appearing at around 100 Mev [28|. If the fraction of the emitted energy which goes into photons is e^. 



e^ripB/f (A/, 



then the density of the radiation at this energy is expected to be VL^ 
observed 7 -ray background density around 100 Mev is fi^ ~ 10^^, we get flpsH < 10" 
With the use of all these parameters, Eqn (27) leads to an upper bound 

I3{M^) < 0.746 X lO"'^ 



Since ~ 0.1 |29f and the 



(28) 



Here for the presently evaporating PBHs the upper limit is much greater than previous results j23l . l30l . |3l[ obtained 
by assuming Einstein's theory or Brans Dicke theory as the theory of gravity. But from our calculation, we found that 
the formation time of presently evaporating PBHs is nearly same in all theories. This higher upper bound implies 
that in LQC a much larger number of PBHs would form in early Universe in comparison with standard cosmology 
and scalar-tensor theories. 



VII. CONCLUSION 



We have studied PBH evolution in loop quantum gravity. We have estimated the cosmic scale factor a(t) and 
energy density p (t) of the universe for both radiation dominated era and matter dominated era. Both expressions for 
a(t) and p{t) are different from those in standard cosmology. In the limit pc 00 these expressions go over to those 
of standard cosmology since standard cosmology envisages nearly zero size for the universe at the time of creation. 
Using these results as inputs we have studied evolution of PBHs using both accretion of radiation and evaporation. 
We find accretion of radiation has no effect on PBH evaporation in the present formalism. From numerical calculation 
it is found that the PBHs created before 1.443 x lO^^s § could evaporate completely in radiation dominated era 
and the accretion efficiency does not affect the evaporation of individual PBHs formed at different times in this era. 
Further, we found that the upper bound on initial PBH mass fraction is much greater than all previous analyses 
but the formation time of presently evaporating PBHs is nearly the same. The greater upperbound implies that a 
large number of PBHs could possibly form in early universe within the context of LQG in comparison with standard 
cosmology and scalar-tensor theories. 
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